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Abstract 

We study the formation and the shape of a liquid meniscus in a wedge with 
opening angle 2 (p which is exposed to a vapor phase. By applying a suitable 
effective interface model, at liquid-vapor coexistence and at a temperature 
T v we find a filling transition at which the height of the meniscus becomes 
macroscopically large while the planar walls of the wedge far away from its 
center remain nonwet up to the wetting transition occurring at T w > T„. 
Depending on the fluid and the substrate potential the filling transition can 
be either continuous or discontinuous. In the latter case it is accompanied by 
a owprefilling line extending into the vapor phase of the bulk phase diagram 
and describing a transition from a small to a large, but finite, meniscus height. 
The filling and the prefilling transitions correspond to nonanalyticities in the 
surface and line contributions to the free energy of the fluid, respectively. 

PACS numbers : 68.45.Gd, 82.65.Dp, 64.70.Fx, 68.35.Md 
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I. Introduction 

Unless much experimental care is provided in growing crystals and cutting them, gener- 
icly on an atomic scale their surfaces exhibit an irregular topography. Numerous theoretical 
[1-10] and experimental [11-14] studies have demonstrated that if such real surfaces are 
exposed to a vapor phase the resulting wetting phenomena [15,16] may differ significantly 
from the corresponding ones on perfectly planar surfaces of the same substrate-fluid systems. 
These studies are focused on the properties of the adsorbed liquid films averaged laterally 
over the statistical irregularities of the substrate topography. 

However, e.g., in the context of structured semiconductor surfaces [17-19], microfluidics 
[20-24], and templates for the selfassembly of small particles [25-27] highly regular nonflat 
lateral surface structures can be produced. But both the theoretical understanding of the 
local wetting phenomena in such structures [28-42] as well as corresponding experimental 
investigations [43,44] are still in their infancies [45]. 

As a paradigm of such structures we consider a single wedgelike cavity with an opening 
angle 2 if and macroscopic extension along its edge; along this y direction the system is taken 
to be translationally invariant (see Fig. 1(a)). The wedge is filled with a simple, volatile fluid 
such that the bulk, i.e., far away from the edge at x = z = 0, is occupied by the vapor phase 
with temperature T and chemical potential ji. The values of /i are chosen such that the 
bulk is close to liquid- vapor coexistence /j,q(T). (Our considerations are equally applicable 
to the two fluid phases of a phase separated binary liquid mixture.) Accordingly, under the 
influence of the substrate potential a liquidlike wetting film will form at the planar surfaces 
leading to a meniscus near the edge. We analyze the height of this meniscus and its shape 
as a function of \x and T. This analysis is based on thermodynamic considerations (Sec. II) 
and on a phenomenological interface model for the local height of the meniscus (Sec. III). 
For a special choice of the effective interfacial potential we are able to calculate analytically 
the shape of the interface in the wedge explicitly (Sec. IV). By considering a suitably chosen 
Landau-Ginzburg- Wilson theory for an order parameter in this geometry (Sec.V) we can 
discuss the range of validity for the effective interface model studied in Sec. III. Our results 
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are summarized in Sec. VI. 



II. Macroscopic description 

In Sees. Ill and IV we shall analyze microscopically the formation of a meniscus of a 
volatile liquid in thermal equilibrium with its vapor phase. Nonethless it will turn out to 
be instructive to consider first in the following subsections the corresponding macroscopic 
description of the problem in constrained equilibrium and to refer to the results obtained 
within the standard theory of capillarity. 

A. Constrained equilibrium and thermodynamics 

We consider a symmetric wedge (i.e., both walls consist of the same type of substrate 
material) with opening angle 2ip < n (see Fig. 1(a)). It is cut off at a macroscopic height H . 
The lower part of the wedge is filled with liquid; the corresponding wall-liquid, vapor-liquid, 
and wall- vapor surface tensions are denoted as a w i, ai g , and a wg , respectively. The vapor and 
liquid phases are taken to be at bulk equilibrium so that they can exchange volume with no 
cost in bulk free-energy. However, in this subsection we consider a constrained equilibrium 
such that the volume V of the liquid phase is prescribed. The liquid meniscus is enclosed by 
the walls at z — \x\ coty? and the liquid-vapor interface is described by f(x) (see Fig. 1(a)). 
Here and in the following we discuss a two-dimensional wedge. Some aspects of the full three- 
dimensional case will be discussed in Subsec.IIB. In the present macroscopic description all 
length scales are much larger than atomic scales so that the free energy F of the system is 
determined by the aforementioned surface tensions only. For a given configuration of the 
interface f(x) one has (see Fig. 1(a)) 

Xi 

F[f] = a lg I dx yjl + P x {x) + 2^—^a wg + 2 a wl (2.1) 
J v smip smif 

-Xi 

df 

with fx — hi Eq.(2.1) the artificial contributions to the free energy generated by the 
cutoff at z = H have been omitted. The liquid-gas interface intersects the walls at x — ± x\ 



3 



with f(±xi) = x\ cot(/9. The surface tensions determine the contact angle 6 via the Young 
equation 

cos@ = a ^ - gtgj (2 2) 

The difference AF[/] = F[/]— 2 rr a W9 /sin y9 of the energies of the wedge with a given amount 
of liquid and of the wedge filled by the gas phase only follows from Eq.(2.1) by subtracting 
the term 2xo(J wg /sm(p. Accordingly the constrained equilibrium profile / minimizes the 
functional 

Xi 



-Xi 



AF*[f] = o lg [ dx \ y/l + P x (x) + — (/(*) - \x\ cotp)\ (2.3) 

J I snap (Jig I 



with the boundary conditions f(±Xi) = Xicoty? and where the Lagrange multiplier A im- 
plements the constraint of constant liquid volume 

Xi 

V = J dx(f(x) - \x\ cotp) . (2.4) 
-x-i 

(In these formulae f(x) is regarded as a single- valued function. In the case that f(x) consists 
of two branches the corresponding separate analysis for each of the two branches leads to 
the same conclusions as the ones described below. Moreover it is necessary to assume that 
+ (f < 7r; otherwise neither the interface nor its branches can be regarded as single valued 
functions.) According to the standard calculus of variations [46] the equilibrium profile 
fulfills the differential equation 

fxx A 



(i + mi °i 9 

with the boundary conditions 

1 ± f x (±Xi) cot ip COS0 



(2.5) 



(2.6) 



(i + m± Xl) y ^ 

Equation (2.6) is equivalent to the statement that the angle between the liquid-gas interface 
and the wall (measured on the liquid side of the interface) is equal to the contact angle 



in the planar case. Thus the equilibrium profile is given by a part of a circle whose radius 
ft — follows from the constant volume constraint (Eq.(2.4)) and intersects the walls at 
the contact angle 6. For 6 < jir — <p, the liquid-gas interface is concave, while for > ^ — (p 
the interface is convex; = \n — <p corresponds to a flat interface. The center of the circle 
is located at (x c , z c ) with 



„ /„ 1 \ COS0 . 

x c = , z c = - R sign + if - - 7T , (2.7) 

V 2 / sini/j 

where the radius R is given by 



R = w(e + v --n+ °»e°»(e + ^)V' . (2 . 8) 

\ 2 sin </? / 

The value of z c is always positive for a concave interface while for the convex one it can be 
either positive or negative, depending on the value of the contact angle 0. The intersection 
of the circle with the wall occurs at 

Xl = -i?sign (q + (p - ^nj cos(0 + ip) . (2.9) 

In the case that the liquid forms a bridge connecting the walls the free energy is given by 

Xi x 2 

F[f\,f 2 ] = a l9 J dxy/l + fcix) + a lg J dx yfl + fl(x) + 
—Xi —x 2 

+ 2 X °- Xl + X2 a wg + 2^^a wl (2.10) 
sm tp sm ip 

where fi(x) and f 2 (x) correspond to the upper and lower liquid-gas interface, respectively, 
while ± Xi and ± x 2 denote the intersections with the walls. The minimization of the con- 
strained free energy with respect to fi and f 2 leads to the conclusion that both interfaces 
are parts of the same circle whose center is given by Eq.(2.7). For each of them the corre- 
sponding contact angle (measured on the liquid side of the interface) has the same value 
as given by the contact angle for a planar substrate. 

The bridge can exist if and only if z c is positive and larger than the radius of the circle, 
i.e., if > \ix + ip. In this case the radius of the bridge has the value 



R = W (20 - 7T - sin2e)~5 . (2.11) 

Thus, depending on the opening angle f and the contact angle 0, the possible shapes of a 
liquid meniscus can be summarized as follows (see Fig. 2): 

> \ 7T + if : bridge 

jty + ip > Q > j7i — ip : single convex interface 

= | n — if : flat interface 

\ n — f > : concave interface . 

For these configurations the difference in free energy between the wedge being filled with a 
liquid of volume V and vapor otherwise and the wedge being filled with vapor only is given 
by 



. / 7T\ /— Z' 71 COS COS(0 + If) \ 2 

AF = sign(e + f--)2a lg W[Q + f-- + — K - , (2.12) 

single interface , 



and 



AF = 4a lg VV (2f + cos 2 © cot fj (2 - vr - sin 2 0) ^ , bridge. (2.13) 

Any configuration with more than one bridge is disfavored by a higher free energy. = 
\tx + ip corresponds to the border case between a bridge and a convex interface for which 
the lowest point of the bridge coincides with the corner of the wedge. 

From Eq.(2.12) one can infer that for > \tx — f, i.e., for a single convex interface or a 
bridge, AF decreases upon decreasing the liquid volume down to the limiting value V — 0. 
On the other hand, for O < \tx — if the free energy decreases upon increasing V, i.e., upon 
filling the whole wedge with liquid. Therefore 

= -7T - f (wetting) (2.14) 
2 

marks a filling transition for the wedge. This case corresponds to a flat interface (see 
Fig. 2) and the corresponding free energy is independent of the volume of the liquid. At 
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fixed temperature, i.e., fixed, a wide wedge is covered only by a microscopically thin 
liquid film which turns into a macroscopic meniscus upon narrowing the wedge beyond the 
critical opening angle ip c = \ix — O. (These considerations can be extended to the case of 
a nonsymmetric wedge. The corresponding results are presented in Appendix A.) If the 
planar walls undergo a wetting transition the contact angle O(T) decreases as a function of 
temperature and vanishes at T = T w , i.e., 0(T = T w ) = 0. 

Thus for a fixed opening angle Eq.(2.14) defines implicitly a filling transition temperature 
T v such that 0(T V ) = ^n—ip. One has T v < T w and T^i^ = T w . This means that the wedge 
undergoes a filling transition at a temperature at which the outer parts of the wedge far 
away from the center still remain nonwet. The same conclusion has been reached previously 
by Hauge from phenomenological considerations based on the Laplace and Young equations 
[34]. Here we obtain this conclusion by analyzing the free energy of the system. 

However, as long as the volume of the liquid drop has a fixed value there is no phase tran- 
sition and the free energy is an analytical function of temperature and of the wedge opening 
angle. Close to (p = \tx — it can be expanded into a series of powers of {tp — (^n — ©))• 
The dominant term of this expansion is 



A F =■ 



2 a lg yV cot p [p - Qvr - e)) . (2.15) 



The aforementioned qualitative change of the interfacial shape and the change of the behavior 
of the free energy upon increasing V are traces of the filling transition which comes into play 
when the fixed- volume constraint is removed and the system is at liquid-gas bulk coexistence 
of a grand canonical ensemble. The minimization of the functional AF[f] instead of the 
functional AF*[f] leads to the following equation for the equilibrium interfacial shape: 

f 

_ § = . (2.16) 

(1 + PJ 

Thus the interface is flat and horizontal. (In the general case of a nonsymmetric wedge the 
interface still remains flat but not necessarily horizontal.) In the case of a nonconstrained 
system the surface free-energy is given by 



F(H) = — [H a wg (T) + H(cosQ(T v ) - cos &(T))a lg (T)} , (2.17) 
cos 

where H is the height of the flat and horizontal interface measured from the edge of the 
wedge. For T < T v the above surface contribution is minimized by H = 0, while for 
T > T v this contribution is minimized by the largest possible value for H. The latter case 
corresponds to a wedge completely filled by the liquid. The relevant surface free energy has 
the following form: 



2H 

cos ip 



(2.18) 

a wg (T) + (cosG(T^) - cos 0(T)) a lg (T) ,T > T v , . 

According to Eq.(2.18) F is continuous at T = T ip but exhibits a break in slope as function 
of temperature. Thus the filling transition in a wedge is — similar to a wetting transition on 
a planar substrate — associated with a singularity in the corresponding surface contribution 
to the free energy. But the relevant structural properties such as, e.g., the shape of the 
microscopically thin liquid film covering the wedge are determined by the line contribution 
to the free energy which is singular at T v , too. This differs from the wetting transition for 
which all relevant structural properties are determined by the surface free energy alone. 

An analogous analysis excludes a bridge as a stable configuration for the unconstrained 
system. Again both interfaces for such a bridge must be flat. If their heights are denoted as 
Hi and H 2 < Hi, respectively, the free energy is given by 

F(H l7 H 2 ) = — [H a wg (T) + H x (cos0(T^) - cos Q(T)) a lg (T) + (2.19) 
cos 

H 2 (cos 0(7;) + cos0(T))^ 9 (T)] . 

Below T v this expression is minimized by H 1 — and H 2 = 0. Above T v the minimum 
corresponds to H 2 = while Hi takes its maximal possible value. 

In the above considerations we have discussed a wetting situation, i.e., the bulk is oc- 
cupied by vapor and the wedgelike substrate prefers the liquid. Our results can be easily 
mapped onto the corresponding drying situation in which the bulk is occupied by liquid and 
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the wedgelike substrate prefers the vapor phase. In this case the filling transition of the 
wedge occurs for 

6 = l - 7T + y? (drying) . (2.20) 

This is in accordance with Eq.(2.14) because liquid and gas are interchanged but the contact 
angle is in both cases taken to be the one of the liquid phase, i.e., in Eq.(2.14) 9 must be 
replaced by n — O. If the planar substrate supports the drying transition, i.e., 0(T — > 
T d ) = it, the wedge will be completely filled by the vapor phase for T > T v where now 
0(T = T v ) = \k + if. Again one has T v < Td and T^i^ = Td. Equation (2.20) states 
that, under drying conditions, for all temperatures corresponding to 0(T) < \tx there is no 
opening angle ip c such that the wedge is filled with vapor for ip < ip c . On the other hand, if 
0(T) > \k there is always a sufficiently small opening angle (fi c (T) = 0(T) — \tx such that 
the wedge is filled up with vapor for ip < ip c . 

This observation may be of relevance for preparing super water-repellent surfaces [47]. 
These are surfaces which exhibit a contact angle n > > | ir if they are very smooth. On 
the other hand, if the surfaces — made of the same material — are prepared such that they 
exhibit a very porous structure one observes an apparent contact angle close to n [47]. In 
the present context this observation is in accordance with a filling of (wedgelike) pores with 
vapor at temperatures at which the smooth planar surface is not yet dry, i.e., < ir. 

Similar considerations hold for wetting conditions, i.e., when the wedge is exposed to a 
vapor phase. Equation (2.14) states that for all temperatures corresponding to 0(T) > | tt 
there is no opening angle ip c such that the wedge is filled with liquid for ip < ip c . On the other 
hand, if O(T) < \ ir there is always a sufficiently small opening angle cp c (T) = | ir — 0(T) 
such that the wedge is filled with liquid for <p < ip c . 

B. Theory of capillarity 

In the previous two subsections we have considered the constrained equilibrium of fluid 
configurations which are translationally invariant along the y direction of the edge of the 
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wedge. Effectively this corresponds to a two-dimensional system. In a three-dimensional 
system, a fixed finite volume V of the liquid must have a finite extension in the y direction 
and thus the shape f(x , y) of the liquid-vapor interface depends also on y. Based on 
the classical Young-Laplace-Gauss capillarity theory the equilibrium liquid-vapor interface 
configurations are surfaces of constant mean curvature meeting the bounding walls with 
contact angle 9. 

There exists a sizable body of literature devoted to the solution of this problem for the 
present wedge geometry (see, e.g., Refs. [48-53] and references therein). The most recent 
account of the present knowledge based on mathematically rigorous results and numerical 
evidence is summarized in Fig. 6 in Ref.[53]. 

For > 1 7r there are two possible liquid configurations depending on the opening angle 
2 (p. For ip > — 1 7r > Oan edge blob forms; this is a part of a sphere such that the 
liquid is in contact with both walls and the edge of the wedge. The shape of the liquid- vapor 
interface is convex. If the opening angle of the wedge is reduced such that p < Q — 1 7r > 
the liquid looses contact with the edge of the wedge and a spherical bridge connecting the 
two walls is formed. The transition between these two distinct configurations occurs at 
(/? = © — 1 7r > 0. This means that a preference of the planar substrate for the vapor 
phase, i.e., > \n implies a prefilling of the wedge with vapor for sufficiently small opening 
angles of the wedge. 

For < 1 7r a tubular bridge between the walls is not possible. For a wide wedge with 
(p > — 0>O one finds an edge blob in contact with both walls and with the edge 
of the wedge. Upon decreasing the opening angle of the wedge at p = \k — > 0a 
transition to edge spreading occurs which persists fory? < \n — Q > 0. Edge spreading by 
the liquid is not possible for > | n. This means that a preference of the planar substrate 
for the vapor phase, i.e., > |7r implies a prefilling of the wedge with vapor for sufficiently 
small opening angles of the wedge. These findings are in full accordance with the free energy 
analysis of the effectively two-dimensional system discussed in Subsec.IIA. 
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III. Description of filling transitions by an interface model 

Within the macroscopic description in Sec. II the type of substrate forming the wedge 
enters only summarily via the surface tensions a wg and a w i. In the actual microscopic pic- 
ture the fluid particles are exposed to the external substrate potential V(x, z) exerted by the 
particles forming the wedge. The resulting full number density distribution p(x, z) of the 
fluid particles can be determined, e.g., either by simulations or by density functional theory. 
However, in view of the considerable numerical challenges by such an approach so far only 
hard sphere fluids confined by hard walls have been studied in such full detail [40,41]. The 
accessible system sizes of the wedges which can be studied within these approaches are also 
severly limited. Moreover, without attractive interactions the hard-body systems do not 
exhibit filling transitions. 



A. Effective interface Hamiltonian 

The study of effective models for the liquid-vapor interface exposed to an effective in- 
terface potential [54,55], which takes into account the competition between the substrate 
potential and the interaction potential between the fluid particles, serves as a reasonable 
compromise between a full-fledged microscopic theory and the purely macroscopic picture. 
Within this effective approach the morphology of the liquid-vapor interface is determined 
by the effective interface Hamiltonian (see Fig. 1(b)) 

V(l(x) simp) — V(loo simp) 



oo 



H[l]= J 



— oo 





( 


' <*ig 





1 + 



df(x) 
dx 



- sj 1 + cot 2 if + 



sin ip 



(3.1) 



Since z = \x\ cot ip is the position of the surfaces of the wedge, the local thickness of the liquid 
film measured vertically equals l(x) = f(x) — \x\ cot ip and thus exhibits a cusp at x = 
with l'(x — > ±0) = =Fcot(/9. Thus the first term in Eq.(3.1) corresponds to the cost in free 
energy due to the increase of the surface area of the liquid- vapor interface relative to its flat 



configuration at \x\ = oo where 



—r— ( x 

dx K 



oo, 



cot ip so that \Jl + cot 2 ((p) 



-J— . The 

sin ip 



second part of the Hamiltonian takes into account the effective interaction V between the 
liquid- vapor interface and the wedge-shaped substrate relative to the configuration l(\x\ — > 
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oo) = loo- V denotes the effective interface potential which we take to be of the same form 
as for a horizontal liquid-vapor interface interacting with the corresponding planar substrate 
[54,55]: 

V(L) = a wl + a lg + u(L) + ApApL (3.2) 

with uj{L — > oo) = 0, A/i = po(T) — p, and Ap = pi — p g where pi and p g are the 
number density of the liquid and vapor, respectively. In Eq.(3.1) V(L) is evaluated in 
a local approximation such that L = I (x) sirup is the local thickness of the wetting film 
normal to the near surface of the substrate; therefore these potential terms are integrated 
with respect to s ^^ ; i-e., along the substrate surface (see Fig.l(b)). Thus Eq.(3.1) neglects 
the additional effective interaction of the liquid-vapor interface with the distant substrate 
surface. This approximation is expected to be valid for a rather open wedge (see below). 
The equilibrium film thickness on a planar substrate minimizes V(L) and yields the actual 
substrate-vapor surface tension a wg = V(loo sin ip) of the planar substrate. Due to the 
subtraction of those terms which correspond to the asymptotic behavior for x — > ±oo the 
Hamiltonian 7i[l] is finite for all configurations compatible with the boundary conditions. 
Therefore TC[l] describes the line contribution to the free energy associated with the linear 
extension L y of the wedge in the y direction. In the present mean-field theory we neglect 
the fluctuations of the interface along the y direction so that / depends on x only and 
the line contribution carries simply a factor L y which has been omitted in Eq.(3.1). We 
emphasize that this consideration of the line contribution to the free energy represents 
the most important improvement over the macroscopic description which considers only 
surface contributions. Within mean-field theory the equilibrium shape / (x) of the meniscus 

minimizes this line contribution: 

d 2 f(x) 

cri g ZMlL _ = V'((f(x) - \x\ cot y?) sin ip) , (3.3) 

, fdf(x)V X ' 

y dx j 

f (0) = , f( X -> ±00) = loo + \x\ COt if . 
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Equation (3.3) is known as the so-called augmented Young equation and has been studied 
by Kagan et al. [39] but not with the view of filling phase transitions; recently these authors 
extended their analysis to study eye-shaped capillaries [56]. Due to the symmetry of the 
system we confine our subsequent analysis to x > 0. In terms of l(x) the boundary conditions 
are l(x — > oo) = simp and l'(x = 0+) = — cot if. Integrating Eq.(3.3) yields 



crig 



^ df(x)^ 
sin if + cos ip ^ x ' 



V V 1 + W J 



r \ 2 



= V(l(x) simp) - V(Z«x> siny?) . (3.4) 



With f(0) = and l = l(x = 0) one obtains from Eq.(3.4) 
AV 1 

= — [V(lo simp) — V(loo simp)] = 1 — simp = v(<p) . (3.5) 

(Tig (Tig 

Equation (3.5) is an implicit algebraic equation for the filling height /o in the center of the 
wedge in terms of the opening angle ip and the wetting properties of the planar system, i.e., 
(Jig, V(L), and thus l^. At liquid-vapor coexistence V(L — > oo) = a w[ + a ig (Eq.(3.2)) with 
^-(L — > oo) = 1 — cos (Eq.(2.2)) so that the condition for the filling transition 1(0) — > oo 
can be expressed in terms of the contact angle 6 of the planar system (Eq.(2.2)): 

cos B(T V ) = siny? or 0(T V ) = ^tt - tp . (3.6) 

is the lowest temperature for which / = oo. Thus our microscopic approach confirms the 
results for T v as predicted by the macroscopic theory in Sec. II. (This conclusion even holds 
if the integrand in Eq.(3.1) is supplemented by a term proportional to the mean curvature of 
the interface.) It is rather satisfactory to see explicitly that a microscopic theory for the line 
contribution to the free energy renders the same value for the filling transition temperature 
T|p as the macroscopic considerations based on the surface free energies. On the other hand 
this is to be expected because at the surface free energy of the wedge is nonanalytic (see 
Eq.(2.18)) so that the line free energy has to follow suit. This is analogous to the fact that 
surface free energies are nonanalytic at bulk transitions. 
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Before we turn to a closer analysis of the filling transition we note that in the special 
case of a wide open wedge, i.e., ip = — e with e << 1 the effective interface Hamiltonian 
in Eq.(3.1) reduces to (see Appendix B) 

oo ( r / ,„/ s \ 2 



«w= /- (^) — v 

— oo ^ L 

V(Z(:r) sin ip) — \^(Zqo sin 93) 



+ 



(3.7) 



sin <£> 





1 


) 


2 n ~ V 



«1 , 



so that the equilibrium profile /(x) = l(x) + |x| coty? is determined by 

oi q sin tp - = V'(l(x) sin tp) , rr > . 
ax" 1 

Upon integration one obtains with — +00) = 



(3.8) 



1 2 / rf/(x) 



= Vf7(x) sin<£>) — V(Zoo sin </?) . 



(3.9) 



Thus for the filling height l in the center of the wedge one finds due to — > ±0) = 
=F cot ip (compare Eq.(3.5) 



AV 



cos 2 <p = v(ip) 



(3.10) 



(Jig 2 

According to the relations |cos 2 y? = (1 — siny?) cos 2 (|(|7r — ip)) and 1 — siny? = 
2sin 2 (|(|7T — ip)) the approximate implicit equation (3.10) for l Q differs from the corre- 
sponding full equation (3.5) only by terms of the order e 4 : v(ip) = v((p)(l + 0(e 2 )) and 
v(ip) = 0(e 2 ). This is also true for the equation for the filling transition temperature T v 
(compare Eq.(3.6)) 



cos 0(T^) = 1 — - cos 2 cp , 



1 /l 



«1 , 



(3.11) 



sin ip + 2 sin ( - ^- n — ip 



Thus we conclude that Eqs.(3.7)-(3.11) are reliable approximations for a wide open wedge. 
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B. Filling height 

In the implicit equation (3.5) and its approximation (Eq.(3.10)) the right hand side v(ip) 
and v(<f), respectively, do not depend on l ] furthermore the left hand side AV/ai g remains 
unchanged upon the open wedge approximation. This facilitates a transparent graphical 
solution for the filling height l as shown in Figs. 3 and 4. As anticipated the filling height l 
is larger than the wetting film thickness on the walls of the wedge far away from the edge 
of the wedge. At gas-liquid coexistence, i.e., A/x = the filling height l diverges smoothly 
for T /" T v in the case of critical wetting and jumps to a macroscopic value in the case 
of first-order wetting of the corresponding planar substrate. At T = T 9 the wetting film 
thickness asymptotically far away from the center of the wedge remains finite. Thus 
within the present model we find that the wedge does indeed undergo a filling transition at 
T v and that the order of the filling transition agrees with the order of the wetting which 
takes place at T w > on the corresponding planar substrate for which 0(T — > T w ) = 0. 

As will be discussed in more detail in the following subsection the filling height l diverges 
along isotherms Ap — > for T > T v . In that case in Eq.(3.5) the asymptotic behavior of 
V(L — > oo) = AfiApL + a wl + (7ig can be used. Together with V(Zoo siny?) = a wg and Eq.(2.2) 
this leads to AV/a Uj — > A/iAp/ sm.tp/ai g + 1 — cos© so that from Eq.(3.5) one obtains 

lo A/x - , T > = : r — . 3.12 

simp Ap A/i 

We note that the form of Eq.(3.12) is valid irrespective of the order of the filling transition 
and irrespective of the range of the molecular forces. The latter enter only indirectly via 
6, (Jig, and Ap. Equation (3.12) is in full agreement with the macroscopic description in 
Ref . [34] and the numerical results above T w in Ref. [33] . 



C. Line tension 

As apparent from Figs. 3 and 4 the implicit equation (Eqs.(3.5) and (3.10)) for the filling 
height Iq yields two or even more solutions. For l < oo, i.e., for T <T {p with defined by 
Eqs.(2.14) and (3.6) or for A\x > 0, the equilibrium solution is that one whose corresponding 
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profile l(x) with 1(0) = l minimizes the line contribution of the free energy (Eqs.(3.1) or 
(3.7)). For l < oo all competing configurations have the same surface free energy. At 
coexistence, i.e., for A/x = 0, lo can become macroscopically large. In that case one has 
to consider both the line and the surface contribution such that if one solution has a lower 
surface energy it wins out irrespective of the behavior of the line tension; the line tensions 
matter only if the surface free energies are equal. (By construction the bulk free energies of 
all configurations are always the same.) Nonetheless, in any case it is interesting to study 
the thermodynamic behavior of line tension. 

From Eq.(2.18) one infers that at coexistence the temperature dependence of the surface 
free energy density is given by a wg + Sa(T) with 5a(T < T^) = and 8a(T > T v ) < due 
to &i g (T) > 0). This implies that for T > T v the filled wedge exhibits a surface free energy 
which is lower than the surface free energy of the unfilled wedge extrapolated to T > T v . 
This holds independently of whether cr wg (T) is an increasing or decreasing function of T. 
Thus we conclude that l = oo for all thermodynamic states (A/x = 0, T > T v ). 

For all other thermodynamic states with A/x > one has Iq < oo and the line contribution 
to free energy (Eq.(3.7)) of the corresponding profile f(x) can be determined explicitly. From 
Eq.(3.9) one has ^ = ^f(^/2/(^/a^smip))^V(l(x) simp) — V^oo simp) for x > (upper 
sign) or x < (lower sign). With ^ = ^ ± cotp for x > (upper sign) or x < (lower 
sign) the insertion into Eq.(3.7) yields for the line tension r\ 

ri = 2V2a lg A (3.13) 

with 

lo 



A= I <U 

c 



V(l sin ip) — V(loo sin p) 



\ a ig \ a ig 



V(l sin p) - V(loo sin p) 



(3.14) 



where lo is a solution of Eq.(3.10) and minimizes V(L). The quantity A, which has 
the dimension of a length and which gives the line tension up to the positive prefac- 



tor 2\^2cri g , is the area between the curves \J (V(l sin ip) — VfJ^, sin p))/ai g , \Jv(p) = 
(V(Iq sin p) — V^oo sin p)))/a lg (see Eq.(3.10)), and I = (see Figs.5 and 6). Note that 
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obviously the curves AV and v((p) intersect at the same position l as V AV and ^Jv(ip). 
From this graphical interpretation one infers immediately that only the solutions l® > 
have to be taken into account; as expected the filling height in the center of the wedge is 
larger than the thickness of the wetting film far outside and these solutions are the ones 
which increase for T /* T v . 

One can infer the filled state for A/i = and T > T v by considering isotherms A/i — > for 
T > T v (see Eq.(3.12)). Figure 7(a) describes the continuous filling of the wedge for T > T v 
and A/i — > in the case of a critical wetting transition. For an underlying first-order 
wetting transition Fig. 7(b) demonstrates that, along an isotherm with T > T v , at A/j, f(T) 
one encounters a thin-thick transition of the filling height which is not accompanied by a 
phase transition in the structure of the wetting film far away from the center of the wedge. 
We call this phase transition a prefilling transition. Once this prefilling transition locus 
has been passed, the filling height l diverges continuously for A/i — > (see Eq.(3.12)). 
This confirms, within the present interface model, the expectation that the thermodynamic 
states (A/i = 0, T > T v ) correspond to a filled wedge. The order of the filling transition at 
coexistence is linked to the order of the underlying wetting transition. 

From Fig. 7(b) one infers that a jump in the thickness 1^ of the wetting film far away 
from the center of the wedge upon crossing the prewetting line enforces a discontinuity 
in the line contribution to free energy associated with a discontinuity of the whole profile 
l(x), however such, that — surprisingly — lo — l(x = 0) happens to change only smoothly. 
On the other hand, upon crossing the prefilling line l changes discontinuously without 
a change in l^. This behavior confirms the general picture that a nonanalyticity at one 
thermodynamic level (bulk, surface, line, ...) induces nonanalyticities at the same locus 
at all subdominant thermodynamic levels and that each thermodynamic level can develop 
new nonanalyticities at loci where all higher thermodynamic levels are strictly analytic: the 
phase boundaries in the bulk free energy are lines of nonanalyticities both for the surface 
and the line contributions, and the prewetting line nonanalyticity of the surface free energy 
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is also the locus of nonanalyticities in the line tension. On the other hand the bulk free 
energy is analytic along the prewetting line and the bulk and the surface free energy are 
analytic along the prefilling line. 

If the wedge fills, on each side of the wedge a three-phase contact line between the 
substrate, liquid, and vapor is forming. In the limit Z — > oo these two contact lines become 
independent and each of them reduces to the structure of a single three-phase contact line 
on a planar substrate. Therefore one expects that in the limit Iq — > oo the line tension given 
by Eqs.(3.13) and (3.14) should reduce to twice the line tension r\ v i an of the corresponding 
single contact line on a planar substrate: 

oo 

r](l — > oo) = 2^2a lg J dl j yV(/ sin ip) - V(loo sin ip) - ^J^cosifj (3.15) 

loo 

where within the present interface model rj p i an is indeed given by [57] 

oo 

Vplan = ^hg J dl yV(l) -V(U " V^} (3.16) 

loo 

with the spreading coefficient S = a wg — o\ g — a w i = —ai g {\ — cos0) equal to — \<Ji g cos 2 (p 
for T = Ty. However, within a full theory one expects that 77(^0 — > 00) = 2r) p i an + rji where 
r/i is the line tension of a wedge filled with liquid, including the liquid phase as boundary 
condition in the bulk limit. But this latter contribution is not contained in the present 
interface model. 

D. Phase diagram 

In order to translate the general features of the line tension discussed in the previous 
subsection into an actual phase diagram one has to specify the functional form of the effective 
potential, i.e., u(L) (see Eq.(3.2)). In the spirit of the square gradient expression (Eq.(3.7)), 
which is applicable for systems with short-ranged forces [58], we choose the generic form 
obtained by Fisher and Jin [59] 
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uj{L) = W exp ^-^j +U (l -CW 2 ^j exp (-y) + AxiApL . (3.17) 

Since at present we are aiming for the generic, possible features of the phase behavior in 
wedges we refrain from studying the specific effects due to an algebraic decay of uj[L — > oo) 
as it is characteristic for actual fluids governed by dispersion forces [15]; this is left to future 
studies. In Eq.(3.17) £ is the correlation length in the bulk of the wetting phase, i.e., the 
liquid phase. 

For C < Eq.(3.17) yields a continuous wetting transition at coexistence A/x = if 
W{T < T w ) < 0, W{T = T w ) = 0, W(T > T w ) > 0, and U(T W ) > 0. In this case the 
phase diagram consists of the line of first-order gas-liquid bulk transitions at /x — Ho(T) and 
of the temperatures T c , T W) and T v as determined by Eq.(3.11) on that line corresponding 
to the bulk, surface, and line nonanalyticities, respectively. For /x ^ fJ>o(T) there are no 
nonanalyticities. 

For C > 0, W(T) > 0, and U > Eq.(3.16) describes a first-order wetting transition 
at T = T w due to a decrease of W{T) for T / T w such that W{T W ) = W W {C, U) > 0. For 
reasons of simplicity we take U, C, and £ as constant and consider a linear temperature 
dependence of W{T) = W w + W (T W - T)/T w , W > 0. The chemical potential difference 
can be expressend in terms of the dimensionless variable (AfiAp£)/W w . Figure 8 illustrates 
the phase diagram predicted by this model. The thick lines indicate the bulk singularities 
at /i — /i , which for simplicity we have taken to be temperature independent, and the 
prewetting line emanating from T w and ending at the prewetting critical point C pw . The 
prewetting line T pw {Ap) (or A/i pw (T)) joins the gas-liquid coexistence curve tangentially 
[60] such that T pw (Ap — > 0) — T w ~ A/iln A/x in accordance with the exponential decay of 
uo(L — > oo). On the present scale of Figs. 8 and 9 this tangential approach is not visible. 
The first-order filling transition occurs at T v < T w such that T v approaches T w for ip — > |7r. 
This infinite jump from a microscopic filling height to a macroscopic height at coexistence is 
reduced to a finite discontinuity off coexistence /x < /xo forming a prefilling line /x p /(T) which 
joins the gas-liquid coexistence curve also tangentially. From our numerical analysis we find 
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fi pf {T \ T v ) - /x ~ a(AT) 2 + 6(AT) 4 + ... with AT = T - T v . The thin-thick jump of 
the filling heigth across the prefilling line diverges ~ (A/x) _1 for /x — > /x (see Eq.(3.12)) and 
vanishes upon approaching the prefilling critical point C p f. For sufficiently narrow wedges 
the prefilling lines are completely below the prewetting line and shorter than the latter. 
Upon increasing the opening angle of the wedge the prefilling line slides into the prewetting 
line thereby breaking up into two pieces: one between and a so-called extraordinary point 
denoted as and another between a second extraordinary point and the prefilling 
critical point C p f. In the limit (p — > these two pieces shrink to zero, such that E^ — > w 
and E^ — * C pw . Figure 9 summarizes the type of the phase transitions in the meniscus 
shape of the liquid in the wedge across the various phase boundaries for a very wide wedge. 
Crossing the pieces of the prefilling line along the paths 1 and 2 leads to a discontinuous 
increase of the filling height of the wedge but does not change the thickness of the wetting 
layer far away from the center of the wedge. A subsequent crossing of the prewetting line 
leads to a discontinuous increase of the thickness of the wetting film and a discontinuous 
change of the meniscus profile l(x) but such that — surprisingly — just the filling height 
l — l(x — 0) increases continuously (paths 3 and 4). Along path 5 both the thickness of the 
wetting film and the filling height change discontinuously; that happens only upon cross- 
ing the prewetting line. According to Fig.9 for suitable model systems along an isotherm 
A/x — > with a temperature T v < T < T w one observes a reentrant prefilling of the wedge. 
For large undersaturations the filling height is very small. It increases until the prefilling 
line connecting C p f and E^ is reached. Crossing it leads to a discontinuous decrease of 
the filling height upon decreasing A/x. A further decrease of A/x leads again to an increase 
of the filling height which jumps to an even larger value upon crossing the prefilling line 
connecting the filling point / and the extraordinary point E^\ Once this prefilling line has 
been passed the filling height diverges continuously for A/x — > (see Eq.(3.12)). 



E. Fluctuations 

Upon crossing the prefilling line the filling height Iq jumps from a value Iq to a larger 
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value Iq so that a volume AV proportional to ^((Iq) 2 — (Iq ) 2 )L tamp is transformed from 
gas into liquid; L is the linear extension of the wedge in y direction. (The above expression 
corresponds to flat menisci at height Iq and Iq , respectively.) Thus effectively ApAV par- 
ticles participate in this phase transition. Therefore in the thermodynamic limit L — > oo 
only at coexistence, i.e., for the filling transition with Iq = oo, this phase transition cor- 
responds to a true two-dimensional system which can indeed support a phase transition at 
T > 0. However, along the prefilling line Iq is finite so that in that case the system is 
quasi-one-dimensional and cannot undergo a true phase transition for realistic interaction 
potentials. Therefore we conclude that the prefilling lines, as predicted by the mean-field 
theory described in the previous subsections, are wiped out by fluctuations in space dimen- 
sions d — 3. Only the filling transitions at coexistence, which can be either continuous or 
first-order, are true phase transitions. These conclusions are in accordance with considering 
the critical points C pw of the prewetting line and C p f of the prefilling line (see Figs. 8 and 
9). Whereas C pw belongs to the Ising universality class in d = 2, C p f would belong to the 
Ising universality class in d — 1 and thus cannot exist. 

One can identify the type of fluctuations which wipe out the prefilling line. If it would 
exist, by imposing suitable boundary conditions at the two ends of the groove for thermo- 
dynamic states at the prefilling line one could generate a stable interface perpendicular to 
the y direction of the wedge which smoothly interpolates between a portion of the wedge 
filled up to Iq and another portion filled up to Iq. However, in this quasi-one-dimensional 
system the fluctuations of the filling height Iq along the y direction, which are not captured 
by mean-field theory, are so strong that for this thermodynamic state at the presumed pre- 
filling line the interface configuration in the wedge breaks up into many domains with Iq 
and Iq, respectively, whose positions fluctuate strongly. 

However, close to liquid-vapor coexistence Iq is very large so that overturning a Iq- 
domain into a Iq -domain and vice versa becomes increasingly improbable. Therefore at the 
prefilling line as function of T or p the filling height Iq(T,p) changes rapidly but smoothly 
between Iq and Iq such that along the prefilling line for p — > p , i.e., Iq — > oo this crossover 
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between Iq and Iq becomes steeper and is confined to a vanishingly narrow region around 
the prefilling line so that the true phase transition at the filling transition at coexistence is 
restored for A/x — > 0. 

The width of this smooth transition region of l (T,fi) at the prefilling line can be es- 
timated on the basis of the finite-size scaling theory for first-order phase transitions [61]. 
These results can be adapted to the present problem following the line of arguments in Sec. 4 
in Ref.[62] where the corresponding smearing out of the prewetting line on a cylindrical 
substrate of radius r has been analyzed. Up to a pre-exponential factor, which depends 
inter alia on the details of the effective interface potential, the temperature range 5 T within 
which Z (T, fx) crosses over smoothly from to Iq is given by 

where R is a numerical factor of order unity and £ is the energy required for the formation of 
a domain wall between a Iq -domain and a Iq -domain. As a crude estimation we approximate 
£ by 2aig (/q ) 2 tan(/9. (Here the effective width 21q tamp of the wedge replaces the cylinder 
radius r in Ref.[62]; as in Ref.[62] the line tension £; introduced there is appoximated by 
IqCFi 9 .) Since is of order unity [62], where £ is the bulk correlation length, we finally 

arive at the estimate 

S lf - ~ exp (-k (f^j tan^j . (3.19) 

As soon as Iq becomes significantly larger than £ the temperature region 5 T for the smooth 
crossover is vanishingly small. Since Iq diverges as (A/i) -1 for A/j — > along the prefilling 
line we conclude that close to liquid-vapor coexistence the difference between a true first- 
order thin-thick transition for l cannot be experimentally distinguished from the actual 
smooth but very steep crossover. In this sense the prefilling line as obtained by mean-field 
theory remains an experimentally accessible line of (quasi-)nonanalyticities. Only close to 
C p f this smearing out of the prefilling line becomes effective. There, in Eq.(3.19) Iq must 
be replaced by Iq — Iq. 
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IV. Meniscus shape 

So far we have discussed the configuration in the wedge only in terms of its key char- 
acteristic feature, i.e., the filling height Iq (see Subsec.IIIB). The more detailed information 
about the full meniscus shape requires to solve Eq.(3.3) or its approximate version given by 
Eqs.(3.8) and (3.9). Whereas the former typically requires a numerical solution, the implicit 
solution l(x) of the latter reduces to an integration: 

[T~ l{x) di 

- \ o a ig sin V? / / = = = x , x > , (4.1) 

V2 / JV(1 simp) - ^(/oosinv?) 

where Iq is the equilibrium solution of Eq.(3.10) and l(—x) = l(x). (Here and in the following 
we drop the overbar which indicates the minimum of Eq.(3.7).) For a given effective interface 
potential V(L) Eq.(4.1) can readily be solved numerically. 

For the model given in Eq.(3.17) it turns out that one can obtain explicit solutions for 
C = and A/i = describing a critical wetting transition of the corresponding planar 
substrate, i.e., W = W (T w - T) /T w where here W < and U > 0. Within this model 
one has 



u = J-J^) = J-JMJ-k) , t = ^l, (4 . 2) 

smip \ I Wq 1 1 J smip \\l Oig cos ip t J l w 



and 



sin ip t — t v ' 

which diverges for t \t tp where 



k = loo + ~ r~ In — ^— , (4.3) 



T w \W \ 



t v = = V ,„; , cosv?. (4.4) 



The contact angle is given by 



a; (Zoo siny?) 
cosB = 1 H . 



(4.5) 



which leads to (see Eq.(3.10)) 
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1 — COS 



11 

2 U, 



cos 2 If 



v(<p) 



(4.6) 



which is in accordance with Eq.(3.11) so that 
"1 



(4.7) 



For x > the profile is determined by (see Eq.(3.9)) 

l(x) 



. diix) fru (w 

soup— — = J — + exp 

ax y aig \2U 



sin (p 



(4.8) 



dUx 

with the boundary condition 



be written as 



l(x) = loo + 



In 



sin (/? 



x 



1 + 



=0+ 



— cot ip. The explicit solution of Eq.(4.8) can 



cxp 



t x 

hi 



cosp> 



, x > 



(4.9) 



V ^ " 1 

Equation (4.9) is in accordance with all expected limiting behaviors: 1(0) = l as given in 
Eq.(4.3); t v — > for ip — > |7r (Eq.(4.4)) so that for x and t fixed l(x) — > ^ which itself 
reduces to the planar value l& = £ In ^ y^^- ^ (Eq.(4.2)); and for large x the film thickness 
Z(x) approaches its asymptote exponentially from above: 



l(x — > OO) = /qo + — 



£ 1 ~t^t cos V 9 
i e r ^ s 



sm</7 ^ _ i 



(4.10) 



provided T is not too close to T v , i.e., ^ 1 >> exp ^— ^— ^cosip^j. For any fixed value 

of x the profile diverges for t — > ^ as (see Eq.(4.3)) 



/(x , i -> O = / ^— § cos 

The maximum curvature, i.e., V"(xq) = occurs at 



(4.11) 



x = ± 



cos 



In 



' i ^ 



-L - 1 



= ± l 



t - t u 



(k - loo) tany? + 0(t - t v ) 



(4.12) 
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Thus for t — > t v the position of the maximum curvature is given by the intersection of the 
asymptote f(x) = + \x\ cot<p and the horizontal z = 1$. Figure 10 illustrates the change 
of the meniscus shape upon approaching the filling transition temperature t v . 
The excess coverage T (see Fig. 10) associated with the meniscus is given by 



oo 



T = 2Ap J dx[l(x) - /J = 




sin 2 ip t 



( 



t 



(4.13) 



y 

with [63] I(y) = J dxx~ l ln(l + x) (see Fig. 11). Since I(y —> oo) = ^n 2 + |ln 2 y one finds 



o 



that upon approaching the filling transition the excess coverage diverges as 

For fixed temperature r — > vanishes as \k — (p (see Eq.(4.4) and I(y — >• 0) = y). 

The line tension associated with the meniscus shape is given by Eqs.(3.13), (3.14), (3.2), 
(3.17), and (4.2) - (4.4). One finds 

' t \ , / t u 



l+ --l ml-^ 



cot <p . (4.15) 



The line tension is negative and approaches its minimal value at t — with a logarithmic 
singularity (see Fig. 12) : 

V (t -> f„) - v (t = t v ) ~ - ^ - l) in (7; " x ) • ( 4 - 16 ) 

For </? — > |7T (i.e., i/^ — > oo) the line tension vanishes as 

77 ^ -> oo^j = -Zo Xg (^-pj coty? . (4.17) 

Thus for fixed temperature rj (ip — > |7r) vanishes ~ Qff — (see Eq.(4.4)). 

Finally it should be pointed out that in terms of the variable t/t v both the excess cover- 
age T and the line tension 77 (see Eq.(4.15) and Fig. 12) can be expressed by scaling functions 
which are independent of ip and the model parameters ai g , Wo, and U . It will be interesting 
to see to which extent this feature is established by more realistic models. 
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V. Derivation of the effective interface Hamiltonian 



The results presented in Sees. Ill and IV are based on the effective interface Hamiltonian 
Ti[f] given by Eq.(3.1). Although this expression for H[f] is rather plausible the status of 
this equation is that of a phenomenological ansatz. In this section we describe the derivation 
of Eq.(3.1) from a more basic model in order to justify our choice of T~L[f] and to gain insight 
into the limitations of this form of H[f]- 

For an actual fluid confined to a wedge the appropriate approach is to start from a 
density functional theory for the inhomogenous number density distribution incorporating 
the full substrate potential of the wedge and the dispersion forces between the fluid particles. 
However, here we refrain from this very demanding task of bridging the gap between a fully 
microscopic description on an atomic scale and Ti[f]. Instead, by aiming at conceptual 
insight we focus on the less ambitious goal to derive TC[f] from a suitably chosen Landau- 
Ginzburg- Wilson theory for an order parameter m which corresponds to the deviation of 
the mean local number density p(r) of the fluid particles from the average bulk densities, 
i.e., p(r) = \ {pi + p g ) + m(r) such that m < corresponds to a gaslike configuration and 
m > to a liquidlike configuration. The values m < (> 0) are favored by the conjugate 
field h < (> 0). This bulk field is proportional to Ap = po(T) — p such that for h > 
there is the liquid phase whereas for h < one has the gas phase. For the present wedge 
geometry (see Fig.l) the natural form of the Landau- Ginzburg - Wilson Hamiltonian is [64] 



Here we have already assumed translational invariance in the y direction. $(m) = &o(m) — 
mh denotes the bulk free energy density. The bulk free energy density at coexistence h = 0, 
i.e., <3>o exhibits two equally deep minima located at m a0 < and > corresponding to 
the gas (a) and liquid (f3) densities, respectively. With the bulk phase being vapor we have 
either h = CT at coexistence or h < off coexistence. In the present context we employ the 




(5.1) 
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so-called double-parabola model which is the simplest form of &o(m) which allows one to 
obtain analytical results [59,65-67]: 

\K£- 2 (m - m a0 ) 2 , m < 
$ (m) = \ (5.2) 

\K^ 2 (m - nifio) 2 , m > . 
£ 7 denotes the correlation length in the bulk phase 7 = a , (5. The parameters are chosen 
such that <&o(m) is continuous at m — 0, i.e., ^~ 1 m Q , = —^mpo. The term — mh in $(m) 
can be absorbed into the double-parabola form as follows: 

i^C 2 (™ - m ah ) 2 + A a , m < 
$(m) = J (5.3) 

\Kif[m - m ph ) 2 + A p , m > , 

where m 7 /j = m 70 + h^ 2 /K and A 7 = —m y0 h — h 2 ^ 2 /{2K) is the equilibrium bulk free 

energy density of the 7 phase. Off coexistence, i.e., for h < the a phase is the stable 

phase with m a h = m a0 + h^ 2 /K as the equilibrium value of the order parameter. The 

requirement of the occurrence of the metastable liquid phase, i.e., > leads to the 

constraint h > —K^ 2 mp . 

For the surface contribution in Eq.(5.1) we choose the common [64] expression 

$i(toi) = -/iimi - ^9 mi (5.4) 

where m,i(x) = m{x , z = \x\ cot</?). The surface field hi is taken to be positive so that it 
has the opposite sign as h. It favors the liquid phase (5 near the wall whereas the bulk field 
h favors the vapor phase a. The surface enchancement parameter g is taken to be negative. 
In Eq.(5.1) the prefactor of the surface term takes into account the increase of the 

surface area as compared to its horizontal projection. In general the LGW-Hamiltonian can 
additionally contain a line contribution corresponding to the edge of the wedge as well as a 
modified surface term including lateral derivatives of the order parameter. Since it turns out 
that these contributions and modifications do not influence the filling transition temperature 
of the wedge we do not consider them here. 
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Following the approach of Fisher and Jin [59,66-68] we employ the so-called crossing 
criterion in order to construct the effective interface Hamiltonian corresponding to the type 
of interface configurations as depicted in Fig.l(b). Order parameter configurations m(x , z), 
which are monotonic and compatible with the boundary conditions to be positive near the 
walls and negative in the bulk, exhibit a line z = f(x) along which m(x,z) is zero. The 
crossing criterion identifies this line as the position of the interface: 

m[x , z = f(x)) = . (5.5) 

The construction scheme for obtaining the effective interface Hamiltonian is to minimize 
Eq.(5.1) under the constraint of the boundary conditions and under the constraint given 
by Eq.(5.5) for a prescribed configuration f(x). The resulting order parameter m(x, z; [/]), 
which corresponds to this minimum and thus is a functional of f(x), is inserted into Eq.(5.1) 
and yields as the surface contribution to Hlgw the effective interface Hamiltonian H[f]. The 
equilibrium shape f(x) of the meniscus minimizes Ti[f]. (In the following we drop overbars.) 
^f|l = yields (see E q .(5.1)) 

( d 2 d 2 \ 

K [ dx 2 + dz 2 J m ^ X ' ^ = $ '( m ^ x ' z ^ ' ( 5 - 6 ) 



and the boundary condition at the substrate 



K f_- 8 gn(x) cot* ^ 



z=\x\cot(p ~ 



hi + gmp 



sin ip 

supplemented by the bulk boundary condition 



(5.7) 



2 = |x| COt ip 



lim m a (x , z) = m ah . (5.8) 

The indices a, (3, and 7 = a, (5 indice that in the special case of the double-parabola model, 
which will be considered henceforth, either the upper or lower part of Eq.(5.3) must be used 
depending on the sign of m; a corresponds to m < and (5 corresponds to m > 0. In 
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Eq.(5.7) we have used the index (5 assuming that m is positive near the wall; otherwise the 
concept of an interface description would not be applicable. 

Following Refs.[66] and [67] the solution of Eqs.(5.5)-(5.8) is expanded into a series of 
terms which contribute increasing orders of derivatives of l(x) = f(x) — \x\ cot ip : 

m 1 = m 77r + m 7 i + m 72 + ... (5.9) 

where the index n indicates the planar limit 2ip = n. For the individual terms one obtains 
within the double-parabola model (Eq.(5.3)) 

m 77r = m 7 o + F j0 (z — \x\ cot <p, l(x)) , (5.10) 

(II 1 / \ \ ^ *^ ^ / \ 

z — \x\ cot (p, l{x)) — — , (5-11) 

and 

( dl(x)\ 2 d 2 l(x) 
m 7 2 = F 7 2i (z - \x\ cot ip, l(x)) I \ + F 7 22 {z - \x\ cot p, l(x)) (5.12) 

The functions F 7i and F l2 i depend on the lateral coordinate x only via the vertical distance 
z — \x\ cotp from the walls and via the local film thickness l(x) but not separately. We 
note that due to the cusp nonanalyticity of l(x) at x = the higher order terms in Eq.(5.9) 
contain distributions which are singular at x = (see, e.g., Eq.(5.12)). Since m(x, z) is 
smooth for z ^ f(x) (see Eqs.(5.3) and (5.5)) these singularities have to cancel each other 
in Eq.(5.9). Since in the following we focus on the lowest terms we do not pursue this aspect 
further. 

The function F^ (u,v) satisfies the differential equation 

/ 1 d 2 s ' 



\sin 2 <p du 2 
with the boundary condition 



i- 2 \ F y0 (u,v) = (5.13) 



K^-F po (u,l(x)) 
ou 



= - {hi + gm/30 + gFf3 (u = 0, l(x))) sinp . (5.14) 



u = 



Equations (5.13) and (5.14) together with Eq.(5.8) lead to the first contribution m 7vr in 
Eq.(5.9) which does not depend on derivatives of l(x): 
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m an (x, z; [/]) = m 



( z - fjx) . N 
1 — exp \ sm(/7 



and 



mp^x, z; [/]) = m m + B + exp j 

B_ exp 

where the coefficients B± are given as 

5+ = 



' z - f(x) 



z - f(x) 



sin(pj + 
sin </? ] 



and 



5_ 



i - gx 2 



t + Qm m X 
1 - £X 2 



with 



X = exp 



Z(x) 



r = 



hi + gm^p 



sin </? 



and 



G = 



+ g 
- g 



(5.15) 



(5.16) 



(5.17) 



(5.18) 



(5.19) 
(5.20) 



(5.21) 



The remaining functions F 7l (u,v), F j2 i(u,v), and F j2 2(u,v) satisfy the following equations: 



and 



/ i d 2 

s sin 2 ip du 2 
' 1 d 2 
sin 2 ip du 2 



d 2 



- C 7 2 F 1± (u, v) = 2 coty? 



<9 2 F 7 o(m, v) 



L 2 F y2 i(u, v) = 2coty? 



dudv 

d 2 F 1 i(u,v) d 2 F l0 (u,v) 



du dv 



dv 2 



(5.22) 
(5.23) 



sin tp du 2 



£ 7 2 F l2 2{u, v) = 2 coty? 



8F^i(u , v) dF l0 (u , v) 



(5.24) 



<9w <9t> 

For constructing the effective Hamiltonian up to square-gradient terms only the functions 
F 70 and F 7l are needed. The solution of Eq.(5.22) is proportional to coty? and only the 
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square of F 7 i contributes to the effective Hamiltonian. The terms proportional to F 7i , 
i = 1,21,22, vanish because these terms are linear in deviations from m 77r which itself 
minimizes the flat substrate Hamiltonian; thus the prefactors multiplying these deviations 
vanish. For a wedge with a wide opening angle, i.e., ip close to one has coty? << 1 so 
that in the spirit, which led to Eq.(3.7), this contribution can be neglected here, too. Thus, 
the effective Hamiltonian is obtained by inserting m 77r as given by Eqs.(5.14)-(5.20) into the 



LGW-Hamiltonian (Eq.(5.1)). For a wide wedge only terms up to 



f dl(x) 
\ dx 



are retained 



leading — after subtracting the surface contribution to the free energy — to 



dl(x) 

+ Xocp{l{x) sin <p) — — cos ip 



dx J >■■■••■ ,•/./• 
V(l(x) sin ip) — V^vo sin ip) 



+ 



sin ip 



The coefficients and Xa/3 and the potential V are given by 



T, a p(l(x) simp) 



K 

sin ip 



dz 



dl 



Xap(K x ) sin V?) = - - 



K 



simp 



x\ cot cp 
/(*) 



dz dl 



oo 

+ J dz 

m 

oo 



dm r 



dl 



dz 



dm a7T dm av 



dz dl 



\x\ cot Lp 



and 



(5.25) 



(5.26) 



(5.27) 



sin ip j dz 
\x\ cot ip 



oo 



V(l(x) simp) = 

M _ * (mah) + _JL_ 



+ simp J dz 
+ $i(m /37r (x, z = \x\ coty?; [/])) 



$(?rw) + 



2 sin ip \ dz 



(5.28) 



In Eqs.(5.26)-(5.28) the profiles m a7r andm^ are given by Eqs.(5.15) and (5.16), respectively, 
and m a h = m a o + ht^/K is the bulk value of the order parameter. The differentiation with 
respect to / is a differentiation with respect to l(x) for any fixed value of x. The stiffness 
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coefficients S Q/ g and Xap as well as the potential V turn out to have the same functional 
dependence on the thickness of the wetting film as for the corresponding planar substrate. 
Here they are evaluated at the local normal distance I (x) sirup to the wall (see Fig. 1(b)). 
For large thicknesses of the wetting film both stiffness coefficients £ Q( g and Xa/3 approach 
constant values which are the surface tension a a p in both cases. The term linear in , 
which is known from studies of wetting on corrugated substrates [67], is multiplied by cosy? 
and thus for a wide wedge it can be neglected as compared with the other two contributions 
in Eq.(5.25). Therefore we can conclude that the phenomenological ansatz for the effective 
Hamiltonian in Eq.(3.7) can be derived systematically from a more basic theory in the limit 
of a wide opening angle of the wedge. The derivation shows which kind of terms are left out 
by the ansatz in Eq.(3.7). 

Moreover, the model parameters entering into Eq.(3.7) can be expressed in terms of 
those of the underlying LGW-Hamiltonian. Within the double-parabola model and the 
identifications a = g and (5 = I one finds 



1 

2 

1 - - > _2 



= nK(&ml + C'mlo) , (5-29) 



and 



w/3 = 7^ K ^p T - hi(r + mpo) - -g(r + m po ) , (5.30) 



W exp ( - + U exp 



with 



W = IKC^m^T (5.32) 



and 



U = K^ l {Qml + r 2 ) (5.33) 
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(compare Eqs.(3.2) and (3.17)). Acccording to Eqs.(5.20) and (5.32) W is negative at low 
temperatures, at which > is large, and vanishes at r = for = hi/\g\ which 
is an implicit equation for the wetting transition temperature. For \g\ > K/^p one has 
U(t — 0) > and the model exhibits a critical wetting transition as studied in Sec. IV. 
Thus if in Eq.(5.31) u(L) is truncated after the first two terms we fully recover the model 
analyzed in the previous section. 

Finally we note that in the limit of a wedge with wide opening angle the resulting form 
of H[f] in Eq.(5.24) does not depend on the special choice of $(m) corresponding to the 
double-parabola model, i.e., Eq.(5.3) but holds for a general functional form $(m), which 
describes two-phase coexistence with a critical point. In that case Eq.(5.13) is replaced by 

F l0 (u , v) = <t>'(F j0 (u , v)) sin 2 if . (5.34) 

Moreover the fact, that in Eq.(5.3) the functions S a/ g, Xa/3, and V exhibit the same func- 
tional form as for the planar substrate, holds also for a general expression for $(m). 



VI. Summary 

We have obtained the following main results for the structure of a fluid exposed to a 
substrate forming a wedge with opening angle 2ip (Fig. 1(a)): 

I. A nonvolatile liquid spreads along the edge of the wedge if its contact angle 6 on the 
corresponding planar substrate is less than — ip (see Fig. 2 and Eq.(2.14)). A vapor 
bubble in a liquid spreads along the wedge if > \ir + if (Eq.(2.20)). Theory of 
capillarity tells that 0(T = T v ) — \ir — if marks also the filling transition temperature 
T v of a wedge by a volatile liquid in eqilibriun with its vapor reservoir (Subsec.IIB). 
At liquid-vapor coexistence of the bulk phases the wedge is completely filled by the 
liquid phase for T >T V although T v < T w where T w with 0(T = T w ) = denotes the 
wetting transition temperature of the corresponding planar substrate; T^i^ = T w . 
The filling transition constitutes a nonanaliticity in the surface contribution to the free 
energy of the liquid confined by the wedge (Eq.(2.18)). 
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2. For T < T v or off liquid- vapor coexistence the surfaces of the wedge are covered by a 
thin wetting film which requires a more detailed microscopic description, e.g., by an 
effective interface Hamiltonian H[f] for the shape f(x) of the ensuing meniscus of the 
emerging liquid-vapor interface (see Eqs.(3.1) and (3.7) and Fig.l(b)). The interaction 
of this interface with the substrate is governed by the effective interface potential V 
(Eq.(3.2)). This description allows one to compute the line contribution to the free 
energy which determines the shape of the meniscus. Without specifying the explicit 
functional form of V the dependence of the filling height l (Fig. 1(b)) on temperature 
and deviation A/x from two-phase coexistence can be discussed graphically (Figs. 3 
and 4). The filling transition at T v and A/x = can be continuous or discontinuous; 
the order of the filling transition is the same as the order of the wetting transition 
of the corresponding planar substrate. Quite generally l diverges ~ (A/x) -1 upon 
approaching coexistence, i.e., A/x — > 0, for T > T v (Eq.(3.12)). 

3. By analyzing the line tension (Eqs.(3.13) and (3.14)) graphically (Figs. 5- 7) one finds 
that a first-order filling transition at coexistence, at which l jumps from a microscopic 
value to a macroscopic one, is accompanied by a prefilling line extending into the vapor 
phase region of the bulk phase diagram (Fig. 8). This prefilling line is the locus of 
nonanalyticities in the line tension; there the surface and bulk contributions to the 
free energy are analytic. Upon crossing the prefilling line the filling height undergoes 
a first-order thin-thick transition (Fig.9). The prefilling line joins the line A/x = of 
the bulk coexistence tangentially. For increasing opening angles T v moves towards T w . 
Accordingly the prefilling line slides into the prewetting line and breaks up into two 
pieces (Figs. 8 and 9) giving rise to rich reentrant prefilling transitions. These general 
features are born out explicitly by model calculations based on a specific choice of the 
effective potential (Eq.(3.17)). 

4. The prefilling transition as obtained from mean-field theory is smeared out by fluc- 
tuations of the local filling height along the edge of the wedge (Subsec.IIIE). Instead 
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of the jump the mean filling height changes smoothly between a small value 1$ and a 
large value Iq near the prefilling line. For small undersaturations A/x the larger value 
diverges ~ (Aii) -1 so that in this limit the temperature resolution required for dis- 
tinguishing between the jump and the actual smooth crossover is experimentally not 
accessible (see Eq.(3.19)). Due to ^(Afj, = 0) = oo the filling transition at coexistence 
persists even in the strict sense. Thus for a system, which exhibits a first-order wetting 
transition in planar geometry, the prefilling line in a wedge should be detectable in 
experiments. 

5. Whereas the thermodynamic behaviour of gross features such as the filling height 
Iq and thus the filling transition itself can be obtained on rather general grounds, 
the determination of the actual meniscus shape requires model calculations based 
on explicit choices for the interface effective potential V. For short-ranged forces 
(see Eq.(3.17) with C=0) exhibiting a continuous wetting transition the meniscus 
shape (Eq.(4.9) and Fig. 10), the excess coverage (Fig. 11), and the line tension (Fig. 12) 
can be obtained analytically. In terms of the reduced temperature variable t/t^ = 
(T w — T)/(T W — T^), within this model the excess coverage and the line tension are 
governed by scaling functions (Figs. 11 and 12) which are independent of the opening 
angle (p and of potential parameters. The scaling functions are nonanalytic for T — > T v 
(see Eqs.(4.14) and (4.16)). 

6. In the limit of a wide opening angle tp the effective interface Hamiltonian T~C[f] 
(Eq.(3.7)) can be deduced from a Landau-Ginzburg- Wilson Hamiltonian (Eq.(5.1)). 
This derivation (Sec.V) points towards additional terms in H[f] which appear for 
smaller opening angles and which are not yet included in Eq.(3.7). 

Acknowledgements : It is a pleasure for us to acknowledge helpful discussions and com- 
ments by T. Boigs, R. Evans, and T. Getta. 
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Appendix A: Constrained equilibrium in a nonsymmetric wedge 

The analysis presented in Subsec.IIA can be repeated for a nonsymmetric wedge char- 
acterized by contact angles 61 and 62 on the left and right side of the wedge, respectively 
Also in this case the liquid with constrained volume V forms a single spherical liquid-gas 
interface or a bridge. In each case the interface intersects the sides of the wedge at angles 
which are equal to the corresponding contact angles. 

In the case of a single interface the radius R of the corresponding circle is 

= re 1 + e, + 2g-, _ /cose, - c»e,\' 

1 2 V 2cosy? J 

COS01 C0S(©1 + Ifi) + COS 2 COS(0 2 + 

+ 2~s^ J (Aij 

while the surface free energy difference between the filled and the nonfilled wedge is given 
by 

A F = ±2a l3 Vvl e ' +e > + 2 *-* - f C ° Se ' - COSe ' X cot^ 

2 \ 2 cos V9 J 

cos©i cos(6i + if) + cos ©2 cos(0 2 + y?)l 5 / A9 x 
+ 2 sinv? J ^ ' 

where the upper sign corresponds to a convex and the lower sign to a concave interface. The 
interface becomes flat and AF vanishes at a temperature T v determined implicitly by 

Qi(T v ) + 6 2 (T^) + 2 = 7T . (A3) 

This temperature marks a filling phase transition for a nonsymmetric wedge if the constraint 
of a fixed volume is removed in favor of a grand canonical ensemble. At T = T v both the 
filled and the nonfilled configuration have the same surface free energies, independent of the 
volume of the liquid. 

A bridge configuration in a nonsymmetric wedge is possible as well. The common radius 
R of both interfaces is 

R = v/yjei + 2 - vr - i (sin2ei + sin2e 2 )| * (A4) 
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while the free energy difference between the liquid bridge and the nonfilled wedge is given 
by 



AF = 4a lg V \<p 



COS 2 61 + COS 2 ©2 



COS 0i — COS S 

2 cosy? 



cot (f 



x jei + 0i - 7T - J (sin 2 0i + sin 2 : 



■>}' 



The bridge configuration occurs provided the following condition is fulfilled: 



'cos0i + cos0 2 s 



'cos0i — cos 2 



> f . 



(A5) 



(A6) 



\ 2 simp J \ 2 cosy? 
For ©i = 2 these expressions reduce to those given in Sec. II. For the symmetric wedge the 
interface is concave if + y? < \ix and convex if + ip > \ix. In the nonsymmetric case it is 
possible for an interface to be convex or concave or even form a bridge if 0i + pi < \k and 
simultaneously 2 + y? > \n- 



Appendix B: Expansions for a wide wedge 

For a wide wedge the opening angle is close to n so that p = \-k — e with e < 1. In 
this limit I , /. is Sm aH for a„ x SO that tb. ** pan of th. inland E,(3,, caa b e 
expanded into powers of / 2 . In terms of l x — f x — cot p , x > one has 

Jl + / 2 - ~r~ = ~~ r~ (fl + 2Z a cosy? siny? + Z 2 sin 2 p) ~ 2 - iX 
v sm y? sm y? I v ' J 

= — (-/ 2 sin 2 y? + Z x siny? cosy? + 0(e 4 )l (Bl) 
sin p 12 J 

= i siay, {/J - cotV} + 0(e 4 ) 

Here we have used the fact that cosy? = sine = 0(e) and that \l x \ is largest for x = with 
\l x (x = 0)| = coty? = tane so that l x = 0(e). This leads to Eq.(3.7). From a systematic 
point of view in the last line of Eq.(Bl) the prefactor siny? = cose can be dropped and 
cot 2 y? can be replaced by e 2 . However, it turns out that it is advantageous to keep the full 
form of these terms (see Eqs.(3.9) and (3.10)). 
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Figure captions 

FIG.l.(a) Macroscopic description of a wedge with opening angle 2p formed by identical 
walls whose surfaces are located at z = \x\ cot <p. The shape of the meniscus is described by 
z = f(x) or l(x) = f(x) — \x\ cot ip. The liquid-gas interface intersects the walls at x — ±x\ 
with a contact angle 0. The system is taken to be two-dimensional and its height is cut off 
at z — H . (b) Same as in (a) on a microscopic scale which takes into account that far away 
from the center of the wedge the meniscus reduces to a wetting film of thickness Zoo sin ip 
covering the walls at large \x\; = l(x — > oo). l = l(x — 0) denotes the filling height of 
the liquid in the wedge. 

FIG. 2. Classification of the equilibrium shapes of a nonvolatile liquid drop (shaded area) as 
function of the opening angle 2ip of a symmetric wedge and of the contact angle 6. In all 
cases the shape of the liquid- vapor interface is a part of a circle. The space within the wedge, 
which is not filled by the liquid, is occupied by vapor. Here the wedge is two-dimensional 
resembling the situation in which a three-dimensional system exhibits translational invari- 
ance along the edge of the wedge. 

FIG. 3. Schematic graphical solution for the filling height Iq according to Eqs.(3.5) and 
(3.10), respectively, for the case of a critical wetting transition of the corresponding planar 
substrate, i.e., A/x = 0. In this case AV exhibits a single minimum at I — which moves 
smoothly to infinity for T — > T w and becomes more shallow for increasing temperature. For 
I — > oo the ratio AV ja\ g attains the limiting value 1 — cosO(T) which vanishes for T — > T w . 
By construction AV is positive. The intersection with the straight line v(p) yields the filling 
height l . There are two solutions 1$ and iff but only the solution l > is compatible 
with the boundary condition associated with Eq.(3.9) (see Subsec.IIIC). For T / T v the 
asymptote 1 — cos©(T) approaches v(<p) from above so that Iq diverges continuously for 
T / TL 
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FIG. 4. Same as Fig.3 for the case of a first-order wetting of the corresponding planar sub- 
strate. For T / the difference 1 — cos0(T) reaches v(tp) and l increases smoothly to 
a finite maximum value at T — T v . At T v there is another solution l = oo. As will 
be shown in Subsec.IIIC, l = oo is the thermodynamically stable solution for T > T v . 
Therefore at T v the filling height l jumps from the finite value to infinity. 

FIG. 5. Schematic graphical interpretation of the line tension in the case of critical wetting 
of the corresponding planar substrate below (a) and at (b) the filling transition temperature. 
According to Eqs.(3.13) and (3.14), for l = the line tension rj is proportional to the ver- 
tically hatched area A, which is to be taken negatively, whereas for I = the line tension 
i] is proportional to the horizontally hatched area, which is to be taken positively. Therefore 
has the lower line tension and is thermodynamically stable. For T /* T v increases, 
the minimum of V AV becomes less steep, and moves out to infinity corresponding to the 
filling of the wedge. Even for dispersion forces with V(L — > oo) ~ L~ 2 the area A and the 
line tension 77 remain finite for T /* T v . The second value of l , i.e., l^f 1 does not correspond 
to solution of Eq.(3.4) or (3.9). 

FIG. 6. Same as Fig. 5 for first-order wetting of the corresponding planar substrate. For 
T / both /qo and l increase but stay finite. At T v the filling height l = 00 is also a 
solution, but the corresponding line tension is larger than that corresponding to the indi- 
cated finite solution by an amount given by the horizontally hatched area, which is to be 
taken positively. For T > T v the surface free energy favors the filled wedge so that lo(T) 
undergoes a discontinuous jump from a finite value at T = T~ to a macroscopically large 
value at T+. 

FIG. 7. Schematic illustration of the filling transition of the wedge along isotherms A/x — > 
above the filling temperature T v for critical (a) and first-order wetting (b) of the correspond- 
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ing planar substrate. Off coexistence \J AV/ <Ji g increases for large / as (A/jApsin tp/ai g )^l^ 
whereas yj AV/ ~ai g attains the finite value y/l — cos for Afi = which is less than yv(jf) 
for T >T H> . In the case of critical wetting there is a single filling height Iq which can possibly 
be thermodynamically stable. For A/j, — > the wetting film thickness increases slightly 
but remains finite whereas Iq diverges continuously. The vertically hatched area diverges also 
so that i](Afi — > 0, T > T v ) — > — oo. In the case of first-order wetting \J AV/ o\ g exhibits a 
global minimum at and a local minimum at l^. Upon crossing the prewetting line of the 
corresponding planar substrate turns into the global minimum. Thus (b) corresponds 
to a case T w > T > T v . For large A/j the area \A 3 \ = —A 3 is smaller then A 2 > so that 
A 2 + A 3 > 0. Consequently ijp corresponds to the global minimum of the line tension. For 
A/i — > the area \A 3 \ increases without limit so that there is a critical value A/x p /(T) at 

(2) 

which | A3 1 = A 2 so that for A/i < A/x p /(T) the filling heigth Iq becomes the globally stable 
configuration. Upon lowering A/x further diverges, as well as A x + A 2 + A 3 , leading to the 
filling of the wedge. For reasons of clarity we have ignored the slight changes in \J AV/ ai g 
for I < upon lowering A/i. A/i p /(T) marks a prefilling transition in the wedge. 

FIG. 8. (a) Phase diagram for the filling of a wedge in the case that the corresponding 
planar substrate exhibits a first-order wetting transition at T w . The thick phase bound- 
aries represent the gas-liquid coexistence curve at fj, — fj, , which for reasons of simplicity 
has been taken to be a straight line, and the prewetting line emanating from T w and end- 
ing at a prewetting critical point C pw . The bulk critical point T c is off the scale to the 
right. We use dimensionless quantities /i* — /Iq = — AfiAp£/W w (see the main text) and 
T* = T/T w so that = 1. The temperatures T Vi denote first-order filling transition 
temperatures which move towards T w for increasing values of </?j. From each filling tran- 
sition point a so-called prefilling line emanates ending in a critical prefilling point C p f. 
The prefilling lines join the gas-liquid coexistence line tangentially as a quadratic function 
whereas the corresponding tangential approach of the prewetting line is logarithmic and 
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not visible on the present scale. For ip — > jir the prefilling line touches the prewetting line 
and breaks into two pieces. For ip = <p 3 in (a) the lower piece between the extraordinary 
point and C p f is shown. On this scale the upper piece between T^ 3 and E^ is not 
visible as well as the distinction between T w and T V3 . This is resolved in (b) on a magni- 
fied scale. These phase diagrams have been obtained for the model defined by Eqs.(3.7), 
(3.10), (3.2), and (3.17) using the following potential parameters: C/W^ = 3.504823, and 
U/W w = 0.197338. For the angles pi = 81.40°, p 2 = 83.12°, p 3 = 84.84°, y? 4 = 85.99°, and 
V? 5 = 87.13° with W Vi = W(T Vi ) one obtains W^/W w = 1.008580, W V2 /W W = 1.005532, 
W V3 /W W = 1.003130, W V JW W = 1.001899, and, W V JW W = 1.000971 so that T* - T*. = 
^ (J^f±- — l^j. Putting numbers on the axis requires to choose a value for the ratio 
For ^ = 1 one obtains T* - T* = 0.008580,0.005532,0.003130,0.001899,0.000971, for 

Wo w ipi 11111 

i=l,...,5 and = 1.594132 and fi*(T Cpw ) - fJ,* = -0.038024. For a different value of ^ 
the temperatures are rescaled linearly according to the formula given above. 



FIG. 9. Types of morphology of the wetting film in a wedge for the various phases within 
a schematic phase diagram. The notation is the same as in Fig. 8. The opening angle (p is 
sufficiently large so that the prefilling line is split into two pieces forming the two extraor- 
dinary points E^ and E^ 2 \ Along the thermodynamic paths 1 and 2 the filling height Iq in 
the center of the wedge increases discontinuously upon crossing the pieces of the prefilling 
line but the thickness of the wetting film far away from the center of the wedge does not 
jump. Along the thermodynamic paths 3 and 4 increases discontinuously while Iq grows 
smoothly upon crossing the prewetting line. Along the thermodynamic path 5 both l and 
loo jump only at the prewetting line. 



FIG. 10. Shape f(x) = f(—x) of the meniscus in units of the bulk correlation length £ in 
the liquid phase. f(x > 0) = l(x) +xcotp with l(x) given by Eqs.(4.2) and (4.9). The 
planar substrate undergoes a critical wetting transition at T w . The temperature is raised 
towards the filling transition temperature T v , i.e., t/t v = (T w — T)/(T W — T v ) — > 1 along 
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two-phase coexistence in the bulk. In units of £ and in terms of t/t v the profile is deter- 
mined uniquely by the dimensionless parameter ^2U/ai g which is chosen to be 2 here. The 
diamonds indicate the position of the maximum curvature of f(x). As explained in the main 
text this position attains a constant distance from the wall for t/t^ — > 1. The shaded area 
corresponds to half of the wedge excess coverage T/Ap (see Eq.(4.13)) for t/t^ = 1CT 2 . The 
dotted line indicates the asymptote of f(x) extended to x = 0. At the present scale the 
temperature dependence of the asymptotes, i.e., of is not visible. 

FIG. 11. Inverse of the reduced excess coverage T* = r(sin(2<p))/(4Ap£ 2 ) as function of 
reduced temperature t/t v (Eq.(4.13)). Upon approaching the filling transition temperature 
T* diverges ~ In (t/t^ — 1) (see Eq.(4.14) and the inset). For <p — > |7r, i.e., t/t v — > oo the 
inverse coverage 1/T* diverges quadratically so that Y ~ Y* / sin(2<p) ~ p — \n for ip — > \it, 
i.e, t/t v — > oo. Note that in terms of the variable t/t v the functional form of Y* is indepen- 
dent of <p and the model parameters <T/ 9 , W , and U . 

FIG. 12. Reduced line tension rj* = rj/(2£ai g cot tp) as function of t/t v (Eq.(4.15)). rj* is 
negative and attains its minimum at t — in a cusplike singularity ~ {t/t^ — 1) ln(t/t v — 1) 
(Eq.(4.16)). For ip — > |7r, i.e., t/t v — > oo, the line tension r/* vanishes as —\tjt^ so that 
rj ~ rfcotp ~ (<p — |7r) 2 for <p — > |7r. In terms of the variable t/t v the functional form of 
rf is independent of <p and the model parameters ai g , W , and [/. 
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